#SCRIPT MODELO MEDIADOR
library(plspm)

#Construção do modelo estrutural (elipses)
distpsic = c(0, 0, 0, 0, 0, 0)
competic = c(1, 0, 0, 0, 0, 0)
respsoc = c(1, 0, 0, 0, 0, 0)
autresp = c(1, 0, 0, 0, 0, 0)
intcomp = c(1, 1, 1, 1, 0, 0)
leacons = c(1, 1, 1, 1, 0, 0)


#Gerar figura do modelo da matriz estrutural
modeloEstrutural = rbind(distpsic, competic, respsoc, autresp, intcomp, leacons)
colnames (modeloEstrutural) = rownames (modeloEstrutural)

#Visualizar modelo estrutural
innerplot(modeloEstrutural)

#Construção dos fatores
modeloMedidas = list(1:3, c(10, 12, 13, 14), 5:9,  15:18, 19:23, 24:28) 
reflexivo = c("A", "A", "A", "A", "A", "A")

#Modelo
resultados = plspm(dados, modeloEstrutural, modeloMedidas, modes = reflexivo)
summary(resultados)

#Visualizar variáveis problemáticas
plot(resultados, what = "loadings")

#Coeficientes de caminho
resultados$path_coefs

#Bootstrap
val = plspm(dados, modeloEstrutural, modeloMedidas, modes = reflexivo,
            boot.val = TRUE, br = 1000)

#Bootstrap results
val$boot
